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Abstract 
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The lecture provides an introduction to thermal field theory and its applications 
\Q ', to the physics of the quark-gluon plasma, possibly created in relativistic heavy ion 

collisions. In particular the Hard Thermal Loop resummation technique, providing a 
consistent perturbative description of relativistic, high-temperature plasmas is intro- 
duced. Using this method interesting quantities of the quark-gluon plasma (damping 
\q ', rates, energy loss, photon and dilepton production) are discussed. Furthermore re- 

cent developments on non-equilibrium field theory, which are relevant for high-energy 
heavy ion physics, are presented. 
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1 Introduction and Motivation 



The aim of relativistic heavy ion collision experiments is the discovery of a new state of 
matter, the so-called quark-gluon plasma (QGP). Quarks are substructure particles 
of hadrons. The observed hadron spectrum can be described by 6 different quark 
flavors (up, down, strange, charm, bottom, and top). Baryons contain 3 quarks, e.g. 
the proton 2 up and 1 down quarks. Meson consists of a quark and an antiquark, e.g. 

7Y + =U d. 

The strong interaction between quarks is described within QCD, where the inter- 
action is caused by the exchange of gauge bosons, the so-called gluons, analogously 
to photons in QED mediating the electromagnetic interaction. In contrast to QED 
quarks have three "charges" (and "anticharges" ) of the strong interaction, called 
color. As a consequence QCD is a non-abelian gauge theory based on the SU(3) 
group. Also in contrast to QED the SU(3) group leads to 8 gauge bosons, which 
can interact directly with themselves. The gluon self-interaction causes asymptotic 
freedom, i.e., quarks and gluons interact weakly at small distances or large momen- 
tum transfers. At large distances, however, the potential between partons (quarks 
and gluons) seems to increase linearly, which explains confinement, i.e., the absence 
of free quarks and gluons in nature. In this picture nucleons and other hadrons can 
be regarded as quark bags, containing besides the valence quarks also virtual quarks 
(sea quarks) and gluons, which force the partons to stay inside of the bag. 

A nucleus can be pictured as a dense system of quark bags. If one increases 
the baryon density, e.g. by compressing the nucleus, or adds further hadrons, e.g. 
by increasing the temperature leading to thermal pion production, these bags will 
overlap. Then the quarks and gluons are not restricted to individual bags anymore 
but can move around in the entire system, which is now in the QGP phase. This 
deconfinement transition is expected to occur at a critical baryon density of the order 
of 10 times nuclear density po, where po — 0.125 GeV/fm 3 = 2.2 x 10 17 kg/m 3 . The 
critical temperature is estimated by lattice calculations (see below) to be in the range 
of T c = 150 - 200 MeV = (1.8 - 2.4) x 10 12 K. These considerations give rise to the 
phase diagram shown in Fig.l. In nature the deconfinement phase transition occurred 
during the expansion of the early Universe about 2 microseconds after the Big Bang 
at zero baryon density (almost equal number of quarks and antiquarks). It might be 
realized also in the core of neutron stars at high baryon density. In the laboratory 
the fireball created in relativistic heavy ion collisions could be in the QGP phase for 
a short time period. 

At this point let me make a few remarks: 

1. The QGP is a relativistic plasma, i.e., the thermal velocity of the up-, down- 
quarks, and gluons is relativistic, since their masses are much smaller than the tem- 
perature of the plasma. 

2. The QGP has the astonishing behaviour that it approaches the ideal gas limit 
for large densities and temperatures, because the parton interaction becomes weak 
due to asymptotic freedom. 

3. Are there alternatives to the QGP at high temperature? Could it be that 
the deconfinement transition is never complete even at arbitrary high temperatures? 
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Maybe new non-perturbative configurations like QCD monopoles exist even in the 
high temperature limit? In order to answer this question the first step must be to 
identify the QGP state in relativistic heavy ion collisions. 

There are two possible scenarios for producing a QGP fireball in relativistic heavy 
ion collisions. At relatively small center of mass energies, i.e. \fs <C 100 A-GeV, 
one expects that the Lorentz contracted nuclei stop each other in the collision and a 
hot, highly compressed fireball with a finite baryon or quark density (finite chemical 
potential), according to the surplus of quarks over anti-quarks from the nuclei, is 
created. At much higher energies, however, a/s~100 A-GeV, there is not enough 
time for stopping and the nuclei penetrate each other transparently. However, the 
vacuum in the space-time volume between the nuclei after the collision is highly 
excited leading to a violent production of gluons and quark-antiquark pairs. Hence 
we expect the formation of a hot parton gas at zero baryon density (zero chemical 
potential). The both scenarios are sketched in Fig.2. 

In the table below former, present, and future experiments with relativistic heavy 
ions are listed. Also a rough estimate for the maximum temperature reachable in 
these collisions is given. 

Accelerator Projectile- Target y/s [A-GeV] T maa , [Me V] 

AGS (BNL) Si-Au, Au-Au 4-5 150 

SPS (CERN) S-U, Pb-Pb 17-20 190 

RHIC (BNL) Au-Au 200 230 

LHC (CERN) Pb-Pb 5500 260 

How does a high energy nucleus-nucleus collision proceed in space and time? The 
space-time evolution of the fireball in the ultrarelativistic case, a/s~100 A-GeV, is 
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Fig. 2 



shown in Fig. 3. The x-axis shows the beam direction and the |/-axis the time, with 
t = at the maximum overlap of the nuclei. The produced particles in this dia- 
gram lie only above the light-cone due to causality. The hyperbolas denote curves of 
constant proper time r = \Jt 2 — z 2 . At the beginning the parton gas will be not in 
equilibrium. Only secondary collisions between the partons will lead to equilibration. 
Afterwards the temperature of the fireball will decrease in the course of its expansion 
until the critical temperature has been reached. Then the system will hadronize, 
maybe showing a substantial mixed phase, depending on the order of the phase tran- 
sition. Finally the system will be so dilute that the interactions between the hadrons 
cease (freeze-out). 

In order to speak of the QGP as a thermal system, we need a large volume and 
particle number, and a sufficient life time of the equilibrated system. The following 
rough estimates from simple models support the possibility of a QGP formation in 
heavy ion collisions: maximum volume (U-U) ~ 3000 fm 3 ; parton number ~ 10000; 
pre-equilibrium period ~ 1 fm/c = 3 x 10~ 24 s; lifetime of QGP ~ 5 - 10 fm/c. 

The main problem is the identification of the QGP phase during the early stage of 
the fireball. Hadronic signatures (e.g. •//^-suppression or strangeness enhancement) 
are often strongly affected by final state interactions in the hadronic phase, while 
electromagnetic signals (photons, lepton pairs) are covered up by a huge background 
of hadronic decays. 

In order to understand the properties of a QGP and to make unambiguous pre- 
dictions about signatures for the QGP formation, we need a profound description of 
the QGP. For this purpose we have to use QCD at finite temperature and chemical 
potential. There are two different approaches. 

1. Lattice QCD is a non-perturbative method for solving the QCD equations 
numerically on a 4-dimensional space-time lattice. In this way all temperatures from 
below to above the phase transition are accessible. Static quantities as the critical 
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Fig.3 

temperature, the order of the phase transition, or the equation of state of the QGP 
have been investigated successfully in this way. However, it is not possible (or at least 
very difficult) to address dynamical quantities, such as the most signatures for the 
QGP formation, a finite chemical potential, and non-equilibrium situations. 

2. Perturbative QCD at finite temperature is based on the fact that the tem- 
perature dependent running coupling constant is small at high temperatures due to 
asymptotic freedom, T — > oo =>> a s (T) = g 2 /4:n — > 0. At a typical temperature 
of T = 250 MeV we expect a s = 0.3 - 0.5. This suggest that perturbation theory 
could work at least qualitatively. It corresponds to an expansion in a s , which can be 
performed conveniently by using Feynman diagrams as the ones in Fig. 4 for elastic 
quark-quark scattering. In this way cross sections, life times, production and decay 
rates etc. can be calculated. The advantages compared to lattice calculations are, 
that one can compute static and dynamical quantities at T > and \i > and that 
an extension to non-equilibrium is possible (see chapter 5). The drawbacks are that 
it is reliable only at high temperatures (T 3> T c ) and that one encounters infrared 
singularities, as we will discuss in the following in detail. 

Finally let me note in this motivation of thermal field theory that there are more 
applications besides the QGP in relativistic heavy ion collisions, namely interactions 
of neutrinos and other particles in Supernovae plasmas, a possible quark matter core 
in neutron stars, the origin of the baryon asymmetry in the early Universe, and Bose 
condensates in condensed matter physics. 
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2 Introduction to Thermal Field Theory 



2.1 Green functions at T > 

The aim of this lecture is the introduction to thermal field theory (TFT) for allowing 
perturbative calculation of cross sections, life times, rates, etc. of particles in an 
equilibrated, relativistic medium. Thermal field theory is a combination of all three 
basic branches of modern physics, namely quantum mechanics, theory of relativity, 
and statistical physics. Therefore one could also call it relativistic quantum statistics. 

Our aim is to derive Feynman diagrams and rules at T > (// 7^ 0). The most 
important quantity in perturbative field theory is the 2-point Green function or prop- 
agator. Therefore we have to ask first: How do propagators look like at T > 0? 
(In the non-relativistic case Green functions at finite temperature are discussed e.g. 
in Fetter, Walecka, Quantum Theory of Many- Particle Systems.) Here we want to 
discuss the scalar field theory first, which we will use in the following as a simple toy 
model to study the different techniques in TFT. 

In order to set the notations and for comparing with the vacuum case, we repeat 
some basic facts about quantum field theory at zero temperature. The Feynman 
propagator is defined as 

iA F (x-y) = (0\T{<i>(x)<f>{y)}\0), (1) 
where the scalar field can be expressed by the Fourier transform 

d 3 k 1 
(2^p Wk) 1 ' 2 



(2) 

k =LO k 



Here x = (xo,x), xq = t using natural units h — c — hs — 1. We use the Minkowski 
metric, x 2 = x^x^ = g^ u x u x fl = Xq — x 2 . Four momenta are denoted by K = (k , k), 
k = |k| and the energy of the field is given by uju = Vk 2 + m 2 . The time ordered 
product of two fields is defined as 



T{cf>(x)<P(y)} 



<f){x)(f)(y) x > y 
<j)(y)<j)(x) x < y 



The Fourier coefficients in (0) represent creation a^(k) and destruction operators a(k), 
which create or destroy a boson with momentum k in the state of the system. In 
particular the vacuum state is given by a(k)|0) = for all k. 
Using ([!]) and (0) the Feynman propagator can be written as 

r d 4 K e -i K <^-y) 

A F (x-y)= ——— — — . 3 

J (27r) 4 K z — m z + ie 

From this expression we can derive the Feynman rule for the propagator in momentum 
space. Ap describes the free propagation of a free scalar particle from y to x for 
^0 > 2/o (creation at y, destruction at x) and from x to y for x < y . Integrating 
over k Q in the complex fc -plane we find in the case xq > yo 

a / x . f d 3 k 1 
A F (x -y) = -% J 



-iK-(x—y) 

(2tt)3 2to k " 



(4) 

ko=u)k 
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Now let us turn to T > 0. Vacuum expectation values have now to be replaced 
by quantum statistic expectation values, i.e. 

(A)=Tr(pA), (5) 

where A is an arbitrary quantum operator and p the density operator or matrix. 
Choosing the canonical ensemble it is given by 

P = \ e-W, (6) 

where (5 = 1/T, H is the Hamiltonian of the system with the eigenvalues and eigen- 
states H\n) = E n \n). (For p ^ H is replaced by H — pN with the number operator 
N .) The partition function is Z = Tr (e~@ H ). Then we can write 

(A) = 1 Tr (Ae-**) = i £ (n\A\n) e"^, (7) 

Z Z n 

where the sum goes over all thermally exited states weighted with the Boltzmann 
factor \n) and exp(— (3E n ). 

(It is possible to formulate the Boltzmann factor in a Lorentz invariant way by 
introducing the four velocity u^. exp (—/3u fJi P fl ) with p = E. In the following we will 
consider only the Lorentz frame of heat bath: = (1, 0, 0, 0).) 

Now we will apply (|7]) to the scalar propagator, yielding 

• A T F >°(x - y) = I £ {n\T{4>{x)<j>{y)}\n) . (8) 

Z n 

Using (0) in (§) for x > y we obtain 

■ r "/ , _ „\ - — f d?k' 1 1 \ - c -pE n 



i A T F >0 (x y) z J (27r)3/2 (27r)3/2 (2 ^ )1/2 (2 ^, )1/2 £ 

(n\[a(k)e~ iK - x + a\k)e iK - x ] [a(k') e - iK '- y + a\k')e iK '- y \ \n). (9) 

The multi-boson states are given by acting repeatedly with the creation operator on 
the vacuum state according to 

|n) = MkO,^),...) = II [ = = |0>- (10) 
The states |n) are orthonormalized. In order to evaluate (|9|) we need 



a(kj)|n) = Jn(kj) |ni(ki),n2(k 2 ), ...,rii(kj) - 1, ...}, 



o + (ki)|n) = v /n(k i ) + l|ni(k 1 ),n 2 (k 2 ),...,n i (k J ) + l,...), (11) 
which leads to 

, . t^>o i < ~\- f d k d k 1 1 % 

1 ^ ~ V) = Z 1 (2vr)3/2 (2tt)3/2 (2w fc )V2 (2c^)V2 ^ e 



{[n(k) + l]5 3 (k - kOe-^ ^'^ + n(k)5 3 (k - k')e iK - x - iK ' v } 

- 1 lwfk^ {[n(v) + 1 ^ ( "" )+ " (k ) e " ;M ) 



(12) 
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Now we use 

\ En(k)e* = 1 = n B {to k ) (13) 

z V ex P (pwfc) - l 

where £7„ = Sk^fc^OO (see e.g. i?ez/ ; Fundamentals of Statistical and Thermal 
Physics) and n B (u k ) is the Bose-Einstein distribution. (In the case of fermions, we 
have according to the Pauli exclusion principle n(k)e{0, 1}, from which the Fermi- 
Dirac distribution np(ujk) = l/[exp(/3(u)k — /f)) + 1] follows. Owing to particle number 
conservation, e.g. charge or baryon number conservation, the average number of 
particles is fixed, which is taken into account by introducing the chemical potential 
ji.) Combining flI2|) and ( p~3| ) we find 

i A T F >°(x - y) = j -j^L JL {[1 + n B (u k )} e~ iK ^ + n B (u k ) e iK ^} . (14) 

Indeed for T = 0, i.e. n B = 0, we recover the vacuum result (|J). 

The physical interpretation of this expression is the following: As at zero temper- 
ature the finite temperature propagator describes the propagation of a scalar particle 
from y to x (x > y ). However, besides spontaneous creation at y there is also in- 
duced creation (~ n B ) and absorption (~ n B ) at x due to the presence of the thermal 
particles in the heat bath. 

Next we want to give a 4-dimensional iT-integral representation of A^ >0 from 
which we can derive Feynman rules in momentum space. 



2.2 Imaginary Time Formalism (ITF) 

We start with the following statement: Going to imaginary times t with < r = it < 
13 and summing over discrete energies k Q = 2ninT (Matsubara frequencies) instead 
of integrating, i.e., 

dko \- 

27r n=-oo 

the propagator can be written as 

""'^^/S^'* (15) 



Proof: (ITBT) can be written as 



d 3 k ^ 1 



Now we use the formula 



T £ f(ko = ^inT) = — / dk - [f(k ) + f(-k Q )} 

Z7TZ J-ioo 2 



1 rioo+e 

+ 7T- / dk [f(k ) + f{-k )]n B {k ), 

Z7TZ J-ioo+e 
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Im k„ 




which holds if f(ko) has no poles on the imaginary axis (see problem #2). 
Choosing 



we get 



with 



-t E 



h 



f(ko) 



-1 



k 2 - (,i 2 



-k r 



-ko-T 



h + h + h + h, 



1 1 



2m 2 

- Res 
2 



dko 



-kor 



too 
-k T 



k 2 - i,) 2 
K ^k 



k 2 - i,) 2 



1 e~^ T 



2 2uj k 

where we have closed the contour in the half plane Re ko > (see Fig. 5). Analogously, 
closing the contour for Re k$ < we find I2 = h ■ 
Furthermore we have 

ioo+e 1 p— fa)T 



2ni J- 



U2 



hi - uj\ e?** - 1 



-UJ k T 



where we have used the contour shown in Fig. 6. For J 4 we have to use the same 
contour as in Fig. 6, since (3 > r, yielding 

p+LU k T 

1 4 = n B (uj k ) 



2oo k 



Combining these expressions we find 



• A^>°(x) 



d 3 k 



ikx 



9.1 j h. V 



2u k 



2u 



-I 



d 3 k 1 



(2tt)3 2uo k 



{[1 + n B (uj k )} e~^ 1 e ik ' x + n B {u k ) e ikx } . 
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Replacing k by — k in the second part of the last expression we recover ([14]). □ 

Another motivation for introducing an imaginary time is the following: for r = 
it — > [3 the Boltzmann factor exp(— f3H) assumes the form of the time evolution 
operator exp(— iHt). As a consequence thermal propagators (more general, all Green 
functions) become periodic with j3, since 

(t x = T,Ty = 0) 

(J3 > r > 0) 



r) 



In general 

A£ >0 (t) = A£ >0 (t + n(3) n integer (16) 

holds. This has two consequences 

1. The time r is restricted to the interval [0, /3[, known as Kubo-Martin-Schwinger- 
or KMS-condition. 

2. The Fourier integral over ko at T = goes over to a Fourier series over the 
Matsubara frequencies ko = 2irinT. 

(For fermions we have Sp >0 (r) = (— l) n Sp >0 (r + n(3) due to a minus sign in the 
definition of time ordering corresponding to anti-commuting fields, from which we get 
k = (2n+l)inT.) 

The Feynman rules in the ITF for example for the 4 -theory now read: 
1. The propagator is given by i A^ >0 (K) = ij [K 2 — m 2 ) with k Q = 2mnT. 



A£>°(x-y) = A-> u (x,y,r,0) 



T>0< 



z 

1 

Tr 



T{0(x,r)0(y,O)} 



e-^0(x,r) 0(y,O) 



z 
z 
z 

A^ >0 (x, y ,r,/3). 



0(y,O) e"^0(x,r) 



0(y,O) e-^0(x, 



-/3H 



0(y,/3) 0(x,r) 



T{0(x,r)0(y,/5)} 
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Fig.7 

2. In loop integrals we have to make the replacement / d A K/ (2tt) 4 — > iT J2k 1 d 3 k/ (2ir) 

3. The vertex reads as in vacuum — i 4! g 2 . (In order to compare with gauge 
theories we denote the coupling constant as g 2 .) 

4. Symmetry factors, e.g. 1/2 for tadpole, are the same as in vacuum. 

As the simplest example for a loop diagram we consider the tadpole of the 4 - 
theory shown in Fig.7. According to the above Feynman rules it is given by 



k 



d 3 k 

% 



(2vr) 3 K 2 - m 



2 



d 3 k ~ 1 



- 12g 2 T [ ^4 E 



(2tt) 3 (2mnT) 2 - u 



k 



1 r d 3 k 00 1 



with Q = u k /(27iT). Using (see Gradshteyn, Ryzhik, Tables of Integrals) 

, . ^ 1 2VL ~ 1 
coth(vTfi) = — — H 



we obtain 



vrfi 7T n 2 + fi 2 



°° 1 °° 1 1 7T 



-oo 



e^n - e-^ f2 e 2nQ - 1 

2vr 2 T 



7T 



( 1 + ^tr)^l 1 + « 



leading to 



ft V e 27rn - 1 / w fc 
r d 3 k 1 

n = 6 ^ 2 /7t4-[ 1 + 2 ^W]- ( 17 ) 

Note that only the vacuum part (ne = 0) is ultraviolet divergent, since the dis- 
tribution functions falls off exponentially for large momenta. Hence we can use the 
same renormalization as at T = 0. Using dimensional regularization the tadpole at 
T = vanishes, resulting in 

47r r°° „ k 2 



Il=12 5( 2 — — / dk — n B (uj k ). 
{2tx) a Jo Uk 
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An analytic expression for this integral is only possible for m = 0, for which we get 
rr 6 2 f°° k x=k/T 6 2rp2 [°°, x 



it 2 Jo e k / T — 1 7r 2 7o e x — 1 

Using 

/*oo ^ n 1 ^2 

/ dx—r=(n-l)\an), C(2) = T 

JO 6 — 1 O 

we end up with the simple result 

II = g 2 T 2 . (18) 

The problem is now that for diagrams with more than one propagator, e.g. the one 
shown in Fig. 8, the summation over ko becomes difficult. A convenient way out of this 
problem is the so-called Saclay method which is a mixed representation performing 
the Fourier transformation in time only. (In the following we use A = A^ >0 .) This 
leads to 

A(T,u k ) = -Tj2e- k0T A(K), (19) 

fco 

where the Fourier coefficients are given as 

A(K) = - / dr e k ° T A(r, u k ). (20) 

JO 

Following the proof of (|15D we can perform the sum (|19|) , 

A(r,w fc ) = {[l + n B (u k )}e-^ T + n B (u k )e^ T }. 

Now the Matsubara frequency ko, over which we have to sum, appears in the propa- 
gator (|20|) only in the exponent, which makes the summation simple at the expense 
of introducing another integral over r. For example for the tadpole we get now 



r d 3 k rP 00 



The sum over k reduces to a (^-function 

oo 

T e ko[T - T ' ] =r^ e 2 ™ T(r - r,) = T5(T(r - r')) = 6{r - t'). 
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This formula makes the Saclay method very convenient in the case of two or more 
propagators. In the case of the tadpole it yields 

n = 1292 Iwr> m " k) 
= 6 «7(S^ 1+2 "^ 

which is identical to our result ([17]). 



2.3 Real Time Formalism (RTF) 



Here we do not aim at a formal derivation, but will present only a motivation and 
plausibility arguments. For a more detailed derivation the interested reader is referred 
e.g. to Ref.[2]. 

The starting point is again our basic equation (|14D for the propagator, 



i A(x — y) 



d 3 k 



(2tt) 3 2cu k 



— {[1 + n B {uJ k )\ e~ iK ^ + n B (u k ) e iK ^} 



k =u) k 



Instead of introducing Matsubara frequencies we will give an alternative expression 
for writing (|T4T) clS Si 4-dimensional integral, namely 



i A(x — y) 



d A K 



K 2 — m 2 + ie 



+ 2tt n B (\k \) 5(K 2 -m 2 



-iK(x-y) 



(21) 



The complex ko- integration over the first term gives the T = term (|j), i.e., the 
first term of ( |T4] ) following from setting n B = 0. The second term gives 



d A K 
d 3 k 



2n n B (\k \) 5(K 2 - m 2 ) e- lK{x - y) 
1 



dk n B (\k \) — [S{ko - u k ) + 5(k + u k )} e -**o(*o-«.) e *-(x-y) 



(2tt) 3 

d 3 k n B (u k ) 
(2tt) 3 2u k 

d 3 k 1 



2uo k 

-iuj k (x -y ) 



+ e 



iu> k {x -y ) 



3 ik-(x-y) 



(2tt) 3 2uj k 



n, B (u k ) e- iK ^+n B {u k ) e iK ^ 



(22) 



k =uj k 



which proves the equivalence of flSID and (p!4|). 

The advantages of this representation are the following. 

1. There is no Matsubara sum, but an integration over ko as at T = 0. 

2. We do not need an imaginary time, which is restricted to finite temperatures 
and cannot be generalized to non-equilibrium. Actually, as we will see in chapter 5, 
the RTF can be used as a starting point for non-equilibrium field theory. 

3. The distribution functions appear from the beginning. Hence the T = and 
T > contributions are disentangled. 
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Fig. 9 



The Feynman rule in momentum space in the RTF for the scalar propagator reads 



now: 



K 2 — m? + ie 

Otherwise we have the same rules as at T 
The fermion propagator is given by 



S{P) = (f + M) 



+ 27in B {\k \) 5{K 2 -m 2 
0. 



(23) 



IE 



2nn F (\p \)5(P 2 -M 2 ) 



(24) 



Again we will consider the tadpole in 4 -theory as an example. In the RTF it is 
given as 



n 



12 g 
12 g 
12 g 



d A K 
d 3 k 



K 2 — m 2 + ie 
1 



dk n B (|A; |) 



(2tt) 3 
o , d 3 k n B (uj k ) 



2u k 



-2n n B (\ko\) 5(K 2 - m 2 ) 
[5(k -u k ) +5(k + uj k )] 



0J k 



where we have neglected the vacuum part, which vanishes after renormalization, 
again. This result coincides with the one, (|i~7|), found in the ITF. 

Unfortunately there is a serious problem with the naive RTF, presented above. 
Singularities, so-called pinch singularities, coming from ^-functions are possible in 
diagrams with two or more propagators, as e.g. in Fig.9. The lower loop in this 
diagram is proportional to / d 4 KA 2 (K) ~ / d 4 K [5(K 2 — m 2 )} 2 ~ 5(0) = oo. 

A possible solution of this problem is the so-called "doubling of degrees of free- 
dom". According to the KMS-condition the time of the fields goes from t — to 
t = —i(3. This contour can be deformed in order to include the real time axis by 
going first from t = to t = oo above the real time axis and then back to t — —i/3 
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below the real time axis. This deformation corresponds to two different kind of fields, 
one existing above and one below the real time axis. Then the propagator, which 
contains two fields, is given by a 2 x 2 matrix (for details see e.g. Ref.[2]) 

A(K) = ( ^jf^ ° )-2m5(K 2 -m 2 ) ( "f°\ k , ^ t, ^ 
V K*-m?-i £ J \ V{k ) + n B {k ) n B (k ) 

(25) 

where we have used here and in the following n B (k ) = ti.b(|A;o|). 

The Feynman rules for the propagator read now i Ay(iiC) with i,je{l,2}. The 
vertex is given by % (— l)- 5 g 2 4!. Fields from above and below the real time axis are 
not mixed at a vertex and the vertex of the "type-2-fields" has an additional minus 
sign coming from the anti-time ordering of the fields below the real time axis. 

Again we consider the tadpole diagram as example. Also self energies are now 
given by matrices. The Iln component, given by Fig. 10, reads 

Hu = \i (-1) 9 2 4! j0p 2vr n B (k ) 6(K 2 - m 2 ), 

which agrees with (f[7|). The components n 12 = n 2 i = 0, since all legs of the vertex 
must have the same index. The component II 22 = —Iln because of the minus sign 
from the "type-2-field" vertex. 

It is important to note that the individual components of the propagator are not 
independent: 

A n - A 12 - A 21 + A 22 

1 - 2m 5{K 2 - m 2 ) n B {k ) + 2m 5(K 2 - m 2 ) [9(-k ) + n B {k )} 



K 2 — m 2 + ie 

+2m 5(K 2 - m 2 ) [9(k ) + n B {k )} - — \ - 2m 5(K 2 - m 2 ) n B (k ) 

1 1 + 2m 5(K 2 - m 2 ) 



K 2 — m 2 + ie K 2 — m 2 — ie 
—2ie 

= 2ilm— — + 2ni5(K 2 - m 2 ) = 0, (26) 

K — m 2 + ie 

where we have used for the last step the important relation 

5(K 2 -m 2 ) = --lm \ (27) 

7r K z — m A + ie 
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1 1 



1 1 



Note that in contrast to ITF factors ie are important in RTF. 

Now we consider the dangerous 2-loop diagrams, for which there are 2 con- 
tributions to rin as shown in Fig. 11. The integral for the lower loop contains 
An (K) An (K) - A 12 (K)A 2 \(K). (Note the minus sign from the "type-2-field" ver- 
tex, which implies that the rules of matrix multiplication do not hold.) Explicitly 
this expression can be written as 



' 2ni5(K 2 -m 2 )n B (k 



K 2 — m? + ie 
x (-2m) 5(K 2 - m 2 ) [9(k ) + n B {k )} 



,K 2 — m? + ie) K 2 — m 2 + ie 
+An 2 5 2 (K 2 -m 2 )[n B (k )+n 2 B (k )] 

(21) ( 1 V 1 i_: w' 1.-2 2, 



- (-2m) 5(K 2 - m 2 ) [9(-k ) + n B (k )} 
Ani 5(K 2 - m 2 ) n B (k ) - 4vr 2 5 2 (K 2 - m 2 ) n% 



K 2 — m 2 + ie) K 2 — m 2 + ie 



AniS(K 2 -m 2 )n B (k ) 



2m \K z — m 2 + te K z — m z — ie 



K 2 — m 2 + ie) \K 2 — m 2 + ie K 2 — m 2 — ie 
A 2 (K) + n B (k ) [A 2 (K) - A* 2 (K)} 



1 + IP? —9 ) 2m 6 ( R2 ~ m2 ) n Bih 



where Ao = l/(-^ 2 — m 2 + z'e). One observes that pinch singularities coming from 
products of 5-functions are absent due to a cancellation between the to diagrams in 
Fig. 11. Similar arguments hold for higher order diagrams. 

Now we will introduce a very convenient representation of the RTF, the Keldysh 
representation. It is constructed from linear combinations of the components of the 
RTF Green functions. The new components of the propagator are defined as 



Retarded propagator: A^ = 


Au- 


Ai 2 , 




Advanced propagator: A^ = 


An- 


- A 21 , 


(28) 


Symmetric propagator: As e 


= An 


+ A 22 . 





The component As is often denoted as Ap in the literature. There are only 3 com- 
ponents, which is sufficient because of ([26]). 
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Fig. 12 



The inverse relations read 

A n 
A 12 
A 2 i 
A 22 



- (A S + A A + A R ), 
^(A S + A A -A R ), 

- (As — Aa + An), 

- (A S -A A -A R ). 



(29) 



Explicitly the retarded propagator is given by 
1 



K 2 — m 2 + ie 
1 

K 2 — m 2 + ie 
f (K 2 - m 2 + ie)- 1 
\ (K 2 - m 2 - ie)- 1 
1 

K 2 — m 2 + i sgn(&o)e 



- 2m 5(K 2 - m 2 ) n B (k ) + 2iri 5(K 2 - m 2 ) [6(-k ) + n B (k )} 

e(-k 1 j 



K 2 — m 2 + ie K 2 — m? 



ie 



k >0 
k <0 



(30) 



It has both poles above the real axis, as shown in Fig. 12. Analogously we get for the 
advanced propagator 

Aa = — 5-^ — (31) 



K 2 — 



,2 — i sgn(A; )e' 

which has both poles below the real axis. The symmetric propagator follows as 



1 



K 2 — m 2 + ie K 2 — m? — ie 
-2ni 5(K 2 - m 2 ) [1 + 2n B (k )}. 



- Aixi5(K 2 - m 2 ) n B (k t 



(32) 



From these explicit expressions, we see that only the symmetric propagator con- 
tains a distribution function. This fact allows for a fast calculation of the thermal 
contributions of diagrams. 

Self energies are related to the bare and full propagators (A, A*) via the Dyson- 
Schwinger equation (see below), II = A -1 — A* -1 , from which together with (Effi) 
follows 

nn + n 12 + n 21 + n 22 = o. (33) 
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Analogously to the propagator we define 

u R = n n + n 12 , 
n A = n n + n 21 , 

n s = n n + n 22 , (34) 

where the different sign in Hr^a is a consequence of (|33|). 

At this point we will give two important relations for the propagators and self 
energies in the Keldysh representation, which we will use later on: 

A S (K) = [1 + 2n B (k )] sgn(fco) [A R (K) - A A (K)}, (35) 

U S (K) = [1 + 2n B (k )] sgn(fco) [U R (K) - U A (K)]. (36) 

The proof in the case of the bare propagator is simple (see problem $4) using explicit 
expressions for the different components ( |30D to (|32|) and the relation (p7|). However, 
it can be shown that (^) holds also for the full propagator in equilibrium. Out of 
equilibrium, on the other hand, modifications are necessary, as we will discuss in 
chapter 5. 

Finally we will consider again the scalar tadpole as an example for the use of the 
Keldysh representation 

. 1 n 2 f d*K 1 A 

= 12 g 2 J 0- 4 27rn B (k o )5(K 2 -m 2 ), 

which agrees with ([171). The last step follows from the fact that the T = O-contribution 
vanishes after renormalization and only contains a T > O-contribution. 
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3 Hard Thermal Loop Resummation (HTL) 

The HTL resummation technique has been developed in the late 80 's and in the 
beginning of the 90's by Braaten and Pisarski (Braaten, Pisarski, Nucl. Phys. B337 
(1990) 569) in order to cure serious problems of gauge theories at finite temperature 
using perturbation theory. For if one uses only bare propagators (and vertices), IR 
singularities and gauge dependent results have been encountered. A famous example 
is the damping rate of a plasma wave in the QGP, which turned out to be different 
in different gauges. 

Braaten and Pisarski suggested the following solution. Instead of using bare prop- 
agators (and vertices) effective propagators, constructed by resumming certain dia- 
grams, the so-called HTL self energies, should be adopted. In this way an improved 
perturbation theory has been invented. 



First we have to isolate the diagrams which should be resummed into effective prop- 
agators. The starting point is the separation of scales in the weak coupling limit. In 
a plasma of massless particles, we have the momentum scale T (hard) and gT (soft) 
assuming g <C 1. HTL self energies are now given by 1-loop diagrams, where the 
external momenta are soft and the loop momenta hard. As an important example we 
will consider the photon self energy or polarization tensor in QED shown in Fig. 13. 
At T = it reads, using standard Feynman rules 



where S is the electron propagator. 

For T > we will adopt the RTF. Then we find according to Fig. 14 for the 
retarded self energy 



3.1 HTL Self Energies 




(37) 




Q=K-P 

Fig. 13 
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Fig. 14 



First we will restrict ourselves to the longitudinal component IT^ = Yl° R . Furthermore 
we neglect the bare electron mass, assuming T ^> m e , which might be realized for 
example in a Supernova plasma (T ~ 10 MeV). We write Sij(K) = Tf. Aij(K), where 
Ay follows from A^- by replacing n B by —n F (compare and (|2"4D). 
The trace over the 7-matrices gives 

tr [Y $ -f%\ = 4 [Q" + K M Q y — gT (Q ■ K)\ 
M= = = °4 (2q k -Q-K)= A(q k + q • k), 

from which we obtain 

U L R (P) = -4z e 2 I ^ (g fc + q ■ k) [A U (Q) A u (tf) - A 21 (Q) A 12 (tf) 

The term in the square brackets of this expression within the Keldysh representation 
reads 

[(A S (Q) + A A (Q) + A fl (Q)) (A S (K) + A A (K) + A R (K)) 
-(A S (Q) - A A (Q) + A R (Q)) (A S (K) + A A {K) - £*(#))] 
= ~ [A S (Q) A R (K) + A A {Q) A S {K) + A A (Q) A A (K) + A R (Q) A R {K)\. 

Note that there are no terms ~ A S (Q) A S (K) ~ 5{Q 2 ) S(K 2 ), which lead to singular- 
ities for P = 0. Also terms A R (K) A A (K), which cannot be integrated, since there 
is no possibility to close the fco-contour (pinch singularity), are absent. 

The /co-integration over ~ A A (Q) A A (K) and ~ A R (Q) A R (K) reduces to zero, 
as one can close the contour always in a half plane that does not contain a pole. Also 
these terms do not contribute to the T > 0-part of the self energy. Hence we are left 
with 

Il£(P) = -2i e 2 J 0^ (q k + q • k) [A S (Q) A r (K) + A A (Q) A S (K) 

Replacing K by — Q in the first term and using A R (—Q) = A A (Q) this expression 
can be simplified further on, 

d 4 K 



K(P) = -4i e 2 / ^ (q k + q ■ k) A A (Q) A S (K) 
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So far this expression is exact. Now we will consider the HTL-approximation, i.e., 
P ~ eT and K ~ T. First we discuss the T = 0-part, corresponding to the 1 in 
the square brackets in the above formula. At T = the only scale is the external 
momentum P. Therefore 11^ _0 ~ e 2 P 2 holds. As we will see below this term is of 
higher order 0(e 4 ) for soft momenta P 2 ~ e 2 T 2 compared to the finite temperature 
part. 

The T > O-contribution yields after integrating over fc 



n£(P) 



16tt e 2 



d A K n F (k ) 



(2tt) 4 2k 
(k -p )k + (k - p) • k 



[S(k — k) + 5(k + k)] 



(k - Po) 2 - (k - p) 2 - isgn(k - p )e 



2tt 3 
■i 



k 



(k-p )k + (k-p) k 



(k - po) 2 - (k - p) 2 - isgn(k - p )e 



(k + Po )k + (k-p) k 



(k + po) 2 - (k - p) 2 - isgn(-k - p )e 



(38) 



Assuming the HTL-approximation, \p \, p <C k, we will expand the integrand for 
small po/k and p/k: 



[...] 



2k 2 — pok — p ■ k 



+ 



2k 2 + p k - p • k 



-2kp + 2k-p + P 2 -ie 2kp + 2k • p + P 2 + ie 



k 



+ 



k 



+ 



1 -Po - W) 



+ ^ 



—Po +pr) - ie po+prj + ie 2 -p + pi] - ie 

1 Vn-Wl 1 P 2 1 P 2 



2 po+pi] + ie 2 (—p + prj — ie) 2 2 (p + pi] + ie) 2 



The integration over 77 = p • k/ (p/c) goes from -1 to 1. The sum of the first two terms 
is odd under i] — > — 77. Hence the first two terms vanish after integrating over 77. Then 
we end up with the final result 



n£(P) 



2tt 2 



dkkiiF(k) / d?7 



-1 



-Po - F7 



+ 



Po-PV 



-po +prj - ie po+PV + it 



(~Po + PV - it) 2 (po+PV + ie ) 2 

o 2/1 Po , Po +P + ie\ 

= -3 m 1 - — In — , 

' V 2p po-p + ie I 



(39) 



where we introduced the effective thermal photon "mass" m 2 = e T /9, coming from 
the /c-integration over the distribution function. The lower limit of this integration 
is zero, although we assumed that k ^> p. However, the error introduced in this 
way is of higher order, as can be seen in the following way. Adopting a lower limit 
eT <^ k* <^ T, we may replace npik) by its zero momentum limit 1/2. Then k* is 
the only scale of the soft loop momentum part of the self energy. Hence it is of order 
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e 2 k* z < e 2 T 2 and can be neglected compared to the HTL contribution. Also note 



that II T>0 ~ e 2 > n 



T=0 



e , justifying the neglect of the vacuum part. 



Analogously we find for the transverse part of the self energy 



n£(P) = 



1 

2 
3 



2 pI 

2 1 p 2 



PiPj 

1 - 



IUj(P) (i, J e {1,2,3}) 



Tr_\ Po_ ln Po + P + ie 
Po) 2 P Po~P + ie 



(40) 



At T = there is only one independent component of the polarization tensor due to 
Lorentz and gauge invariance: 



It 



T=0 



p2 



n(p). 



At T > Lorentz invariance is broken, since we have chosen the heat bath as the 
rest frame. But transversality, P^Yi^ = 0, still holds as a consequence of gauge 
invariance. This leads to two independent components, e.g. U L,T (see e.g. Ref.[l]). 
The advanced self energy follows from the retarded one simply by the replacement 



U L / = U^(ie 



L.T i 



-ie). 



(41) 



Finally, the symmetric self energy is given by (see problem # 5) 

Ate 2 



n L s (P) 



nl(P) 



Tip 



poo 

d{p 2 ~ Po) / dk k 2 n F (k) [1 - n F (k)\ 
J o 



T 



= -Qirim - 6(p -p Q ), 



—3iri 



P 
T 

P 



1 - 



0(p 2 -pD- 



(42) 
(43) 



Note that it is of lower order, 0(eT 2 ), than the retarded and advanced self energies 
for soft momenta p ~ eT. Furthermore it is purely imaginary. 
Now a couple of remarks are in order: 

1. Using the ITF one obtains the same results. In this case one has to continue 
analytically the discrete imaginary energy p = 2mnT to real continuous values. 

2. The HTL approximation is equivalent to the high temperature approximation, 
T ^> p, \p \, which has first been studied by Klimov and Weldon in 1982, and to the 
semiclassical approximation discussed already in 1960 by Silin. 

3. In contrast to the tadpole U L,T are momentum dependent. 

4. For pi < p 2 the self energy has an imaginary part 



p +p±ie 

In = In 

Po - p ± ie 



Po + P 



Po-P 



Tin9(p 2 -p 2 ), 



(44) 



describing the collisionless energy transfer from the plasma modes (see below) to the 
thermal particles of the plasma (Landau damping). 
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Fig. 15 



5. In the static limit, po — > 0, the longitudinal self energy reduces to 

ni x (po 0,p) = -3m5, (45) 



which leads to Debye screening of the electric interaction due to the presence of 
charges in the plasma. However, the transverse part reduces to 



ni A (po-o,p) = o, 



(46) 



i.e., there is no static magnetic screening in the plasma. 

6. Considering QCD instead of QED, the gluon self energy is given by Fig. 15. 
Surprisingly fl39|) to ( |43| ) still hold within the HTL approximation if we simply replace 

9 2 T 2 ( N f \ 

(47) 



m~ 



N f 
1 + — 
6 



where Nf is the number of active flavors in the QGP. Since the 1-loop QED polariza- 
tion tensor is gauge invariant, the same holds for the QCD one, which has the same 
form. 

7. The fermion self energy can also be derived within the HTL approximation 
(see Fig. 16). Starting from the most general ansatz for massless fermions at T > 0, 



Sr(-P) = -a(po,p)f- &(Po,Pho 



(48) 



with 



a(po,p) 
b(po,p) 



— [tr(f Efl)-potr(7oE fl ) 



4p2 



P 2 tr( lo Z R )- Po tr(P E R ] 



(49) 



K~T 



P~gT 




Fig. 16 
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Fig. 17 



one finds in the HTL approximation 



tr(fSii) = 4 m 2 ,, 

tr(7 E R ) = 2m F - \n — (50) 

P Po ~ P + i£ 

with the effective fermion mass m 2 F = e 2 T 2 /8 (QED) and m 2 F = g 2 T 2 /6 (QCD). 
Again one can show that the fermion self energy, (BED to fl5"0"|), is gauge independent. 



3.2 Effective propagators and dispersion relations 
3.2.1 4 -theory 

The pole of the bare propagator, A(K) = 1/{K 2 — m 2 ), describes the dispersion 
relation of a non-interacting, scalar particle: K 2 — m 2 = k^ — io\ = 0, i.e. ko = = 
\/k 2 + m 2 . (Here we use the notation A (if) = Ar(K) and omit ie.) 

We can construct an effective propagator by resumming the self energy using the 
Dyson-Schwinger equation as shown in Fig. 17. This diagrammatic equation reads 

tA* = iA + iA(-iIIi)iA* 
A* = A + AnxA + An 1 AH 1 A + ... 

oo 

= A £(11! A)* 

n=0 

A 1 - n - ! A = a- 1 - n x = K 2 - m 2 - Ui ' ^ 



9<1 



The pole of the effective propagator determines the dispersion relation of an inter- 
acting collective mode with the effective mass M = \/m 2 + H±. (Effective masses, 
generated by the interaction with a medium, have been introduced in various prob- 
lems in physics, e.g. the effective mass of electrons in crystals.) The dispersion rela- 
tion of the collective scalar particle is simply given by the one of a massive particle, 
u(k) = \/k 2 + M 2 (see dashed line in Fig. 20). 
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/x/^'N/x = /\/\/\ + /\A K/\#^/\ 

3.2.2 QED 

In gauge theories we have to fix the gauge in order to determine the dispersion rela- 
tions for a gauge boson from its propagator. At finite temperature it is convenient to 
choose the Coulomb gauge V ■ A = 0. Since Lorentz invariance is broken, the choice 
of a non-covariant gauge is no problem. The bare propagator is then given by 

D 00 (K) = D L (K) = ^, Ah = Ao = 0, 

DiAK) = (^rf)w, W = p (52) 

where we used D = Dr,a, n = T1r,a- The longitudinal propagator Dl cannot be as- 
sociated with a real photon but with the Coulomb potential, i.e. electric interactions. 
The transverse propagator Dt, on the other hand, describes 2 transverse massless 
photons and the magnetic interaction. 

The Dyson-Schwinger equation, shown in Fig. 18 gives the effective propagator 

Dl(K) 



(See problem # 6). 

Using the HTL-approximation for Hl,t we find in the static limit: D^(k = 0) = 
l/(/c 2 — m 2 D ) with the Debye mass m 2 D = 3m 2 ,, leading to screening of electric fields 
by the presence of charges in the plasma. For the transverse propagator we get in 
this limit D^(ko = 0) = 1/k 2 , i.e., there is no screening of static magnetic fields. 

The pole of the effective propagator, D* l.t{K)~ 1 = 0, determines the HTL dis- 
persion relations. Due to the complicated momentum dependence of the HTL self 
energies only numerical solutions are possible, which are shown in Fig. 19. Let me 
make the following remarks: 

1. There are two branches. The longitudinal mode ujl is called a plasmon as 
in non-relativistic plasma physics. It corresponds to a collective longitudinal photon 
mode that is absent in vacuum. 

2. For k — > we have cul(0) = ujt{0) = m 7 which is called the plasma frequency. 

3. For k — > oo we find 0JL,r{k — ► oo) = k, i.e., we recover the free dispersion 
relation, since the self energy II ~ eT can be neglected for large momenta k 3> eT. 
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4. Because of ujl,t > k we have ImlT^T = 0, i.e., the HTL dispersion relations 
are undamped (see below). 

5. The gluon dispersion relation follows from the one of a photon by simply 
replacing m 7 by m g . 

6. In the case of massless fermions the HTL resummed fermion propagator leads 
also to two branches as shown in Fig. 20. The uj + branch corresponds to a collec- 
tive electron or quark mode, whereas the u>- branch, which has a negative ratio of 
helicity to chirality, does not exist in vacuum. Interestingly this so-called plasmino 
has a minimum at finite momentum, leading to interesting consequences for possible 
observables of the QGP (see below). 

Finally we will discuss the symmetric HTL resummed propagator. We start from 
(|35|), which also holds for the full propagator in equilibrium, 



Df{K) = [1 + 2n B (k )} sgn(A; ) [D R L (K) - D\ L {K)}. 
Combining ( 41 ) and (|44]) we see that Imll^ = -Imll^ holds. Therefore we have 

Dr ^ K "> = k 2 -ReU L R (K) ±ilmIL%(Ky 
if Imn|j ^ (i.e., for k 2 > k\ in HTL-approximation) . From this we find 

9?Tm T\ L 

Now we introduce the spectral function pi = —ImD^/n. Then we can write 

D* L (K) = -2m [1 + 2n B (k )] sgn(k ) p L (K). (54) 
The spectral function is of the Breit-Wigner form 

1 imnf 



Pl(K) 



vr (k 2 -Ren^(K)) 2 + (ImU^(K)y 



describing quasiparticles with finite width. Note, however, that the quasiparticle 
"mass" Ren and "width" Imll depend on momentum and energy. In the HTL 
approximation the self energy has an imaginary part only below the light cone (k 2 > 
&o). Then the spectral function can be decomposed into a pole contribution and a 
cut contribution 

p L {K) = pl ol \K)+pT(K), 
p L olc (K) = sgn(k Q )5(k 2 -ReUf(K)), 

„cut m _ _1 Im n «W B(k 2 -k 2 ) 

Pl [ ] ~ 7T (k 2 -ReU L R (K)) 2 + (lmIl L R (K)) 2U[ oh 

where one has to use the HTL results for ReH R (K) and Imn^.K'). 
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Fig. 21 



3.3 HTL resummation technique 

As we said before, naive perturbation theory at T > suffers from a problem, namely 
infrared singularities and gauge dependent results. The reason for this undesirable 
behaviour is that the naive perturbative expansion is incomplete at T > 0. Infinitely 
many higher order diagrams with more and more loops can contribute to lower order 
in the coupling constant. These diagrams can be taken into account by resummation. 



3.3.1 Massless 4 -theory 

The lowest order tadpole shown in Fig. 7 has been calculated already above, IT = 
g 2 T 2 = M 2 . To next order of the scalar self energy we find a diagram as in Fig.9. 
This diagram exhibits a logarithmic infrared singularity because of n 2 ~ g 4 J d 4 K/K 4 . 
Naively we expect that the self energy can be written as II = IT + n 2 + ... = 
g 2 T 2 + 0(g 4 ). 

Now we consider the tadpole diagram II* in Fig.21, where we have an effective 
propagator, iA* —ij [K 2 — M 2 ), in the loop instead of a bare one. Since the effective 
propagator follows from a resummation of the bare tadpole, II* is given by Fig. 22, 
i.e., it is a sum of all daisy diagrams. In order to calculate II* we simply have to 
replace the bare mass in IT by the effective mass M and obtain 



IT = 12 g 2 J — [1 + 2n B (Q K )}, Q k = Vk 2 + M*. 

The fc-integral cannot be done analytically, but an expansion for small g is possible 
(see problem # 8): 

W =g 2 T 2 [l--g + 0(g 2 )}. (55) 

7T 

This surprising result deserves two comments: 




Fig. 22 
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Fig. 23 



1. IT is infrared finite, although it contains (infinitely many) infrared divergent 
diagrams. 

2. The correction to 111 is of order g not g 2 , i.e., it is a non-perturbative correction, 
since it is not a power of the perturbative expansion parameter g 2 . 

These observations suggest the following recipe: 

1. Isolate terms ~ g 2 T 2 . Here: III. 

2. Construct the effective propagator by resummation: Here: A*. 

3. Use A*(K) as in perturbation theory, if k and k are soft (~ gT), because 
then all terms in the denominator of the effective propagator are of the same order, 
A* = l/(g 2 T 2 ). If ko or k are hard (~ T), however, it is sufficient to use the bare 
propagator, which is then of order A = 1/K 2 ~ 1/T 2 . 

In n* we have integrated over k and k from to oo. The soft momentum range 
is important for the correction to III . Therefore it is necessary to use A*. 

3.3.2 Gauge theories 

Here we do not want to discuss the HTL resummation technique for gauge theories 
in detail, but will only take over the arguments from the scalar theory. Our findings 
above suggest the following strategy: 

1. step: Isolation of self energies ~ g 2 T 2 

These are the HTL self energies Hl,t and S, which we discussed already in detail. 
In contrast to the scalar theory there is a new aspect: due to Ward identities self 
energies are related to vertices, e.g. 

le [S(P0 - E(P 2 )] = (P 1 + P 2 )^(P l , P 2 ). 

Therefore we have to consider also HTL vertex corrections, as shown in Fig. 23, where 
all internal lines are hard ~ T . 

2. step: Effective propagators and vertices 

Effective propagators are constructed as above making use of the Dyson-Schwinger 
equations. The HTL vertex following simply by adding the HTL correction to the 
bare vertex. Examples are shown in Fig. 24. 
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Fig. 24 

3. step: Effective perturbation theory 

We use effective propagators and vertices if all legs are soft; otherwise bare prop- 
agators and vertices are sufficient. In this way, contributions of the same (or even 
lower) order in g are included, gauge independent results are obtained, and screening 
effects are included leading to an improved IR behaviour. 

Summarizing, a large progress compared to naive perturbation theory has been 
achieved, although there are still problems, as we will see below. 
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4 Applications 



4.1 Muon damping rate 

The simplest application of the HTL resummation technique for gauge theories is 
the interaction or damping rate 7 of a heavy fermion with mass M in a relativistic 
plasma, e.g. muons in an e + -e~-plasma at m e < T < M M . The interaction rate is 
related to the mean free path by A = I/7. 

The damping or interaction rate describes the damping of a particle with the 
time evolution exp(— iut) (plane wave). In general the frequency u has a real and 
an imaginary part: u = Keu + ilmu. Defining 7 = — Imu we have exp(— iut) = 
exp(— iReut) exp(-'-ft). As an example we consider a scalar field. Its dispersion 
relation is given by u 2 — k 2 — H(u, k) = 0, from which we find 

(Re u - i-f) 2 -k 2 - Re IT(Re u - ry, k) - ilm IT(Re u - ij, k) = 0. (56) 

In the case of no overdamping, 7 «C Reu, we get 

Re 2 u - k 2 - ReIT(Re^, k) = 0, (57) 

from which the dispersion relation Re a; = u(k) follows, and 

-2iReu-f = ilmU(Reu,k) } (58) 

which leads to 

1 = ~2u\k) lmU{uj(k) > k) - (59) 
In the case of a fermion with mass M the interaction rate is given by 

r(£) = ~ [1 - n F (E)] tr [(f + M) Im E R (E, p)], (60) 

where E 2 = p 2 + M 2 and is the retarded self energy of the massive fermion. 

In naive perturbation theory the self energy to lowest order is given by Fig. 25, 
which has no imaginary part on-shell, Im £#(£?, p) = 0. (A bold line denotes a 
massive fermion.) For according to cutting rules the imaginary part of £ can be 
related to the matrix element of Fig. 26, which describes the emission or absorption of 
a photon from a bare muon. However, this process is forbidden by energy-momentum 
conservation. 

To next order we have to consider the 2-loop diagram of Fig. 27. Using cutting 
rules it corresponds \ie ± — > iie ± scattering, which leads to an energy loss of a muon in 




Fig. 25 
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the e + -e~-plasma. From naive power counting we expect T{E) ~ a 2 . However, due 
to the exchange of a massless photon the damping rate turns out to be quadratically 
IR divergent. 

Therefore we have to adopt the HTL method, which leads to the diagram of Fig. 28. 
Due to the non-vanishing imaginary part of the photon self energy contained in the 
HTL photon propagator of the exchanged photon this diagram has a finite imaginary 
part and corresponds to scattering via the exchange of a soft photon (virtual Landau 
damping). Note that the HTL photon self energy contains hard (thermal) electrons. 
Furthermore the HTL resummed photon propagator contains Debye screening. Note 
also that we do not need an effective muon-photon vertex or an effective muon prop- 
agator as the momenta P, P' are always hard due to the large muon mass. 

Using standard Feynman rules we find at T = 




(61) 



At finite temperature, using the RTF, we find 
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Defining S(P') = (f + M) A(P') and using 

D; U Y (f + M) Y^Di^if + M) 7 ° + D* T - 7 l (f + M)^ 

and 

tr [(P + M) T °(f ' + M) 7 °] = 4(2p Po - P ■ P' + M 2 ) 
c% - tr [(P + M)Y(f + M) 7 J ] = 8[p 2 - p q - (p ■ q) 2 + p ■ q - M 2 ] 

we find 

tr [(P + M) 5£(P)] = 

4*e 2 J ^{A n (P) [^(Q) (p 2 + p 2 - p q - p ■ q + M 2 ) 
+2 (p 2 - j9 g + P • q - (P • q) 2 - M 2 )] - [(11) - (12)]}. 

Using the Keldysh representation we obtain 

A n D* u - A 12 D* 12 ^ 1 - [A R D* s + A R D* A + A s D* R + A A D* R }. 

The leading contribution in the coupling constant e comes from n R {kQ ~ eT) = 
l/[exp(|/co|/P) — 1] — P/|^o| ~ l/ e 5 which appears only in the first term. Therefore 
to lowest order we have 

tr [(P + M) S^(P)] ( ^4vre 2 / ^ [1 + 2n B (q )} sgn(g ) 

[(pg + P 2 ~ PoQo - P • q + M 2 ) p L (Q) + 2(p 2 - p go + P ■ q - (p ■ q) 2 - M 2 ) p T (Q)] A R (P'). 
The term ImS comes from Im Aij(P') = — 7rsgn(p ) 5(P' 2 — M 2 ). Hence we find 

T(P) = ^ [1 - n F (E)} J ^ [1 + 2n B (g )] sgn(g ) 

[(2P 2 - Pg - p • q) p L (Q) + 2(p 2 - Pg + p • q - (p • q) 2 ) p T {Q)\ 

sgn(P - g ) ^ [5(E - g - E') + S(E - q + E% (62) 

Now we want to make the following approximations: 

1. E ^> T from which Uf{E) ~ follows. 

2. The HTL approximation for the photon propagator: q, go <C T <C E. This 
leads to 

• [l + 2n B (g )]sgn(g ) ~2P/g , 

• sgn(P - g ) = +, 

• the following simplification in the argument of the first 5-function: 

E' = ^p' 2 + M 2 = y/(p - q) 2 + M 2 ~ ^E 2 - 2p • q 
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• 8{E - g + E') = 0, 

• 1/E'~1/E. 

Using these approximations we find 

T(E) = 4vr 2 e 2 T J I [ Pi (Q) + (t; 2 - (v • q 2 ))p T (Q)] % - v • q) 

Introducing the angle rj 

p • q p 

v ■ q = wgr/, = , v = — 

pq E 

we can use the ^-function for integrating over this angle, using 



8(q -V • q) = — 5(r) - — ), 
t>g t>g 



which gives 



Pl(?o, 9) + f^ 2 - j pt(9o, 9) 



, , . , , , , > , , (63) 

Note that we integrate from to oo, although the momentum q of the photon is soft. 
This is possible since the intergrand drops like 1/g 3 , i.e., the integral is dominated 
by soft momenta regardless of the upper limit as long as it is much larger than eT. 
Hence the hard integration range contributes only to higher orders that have been 
neglected anyway. 

From the intergation range of go we read off that go < g. Hence only the cut 
contributions (Landau damping) of the spectral functions pl,t are needed. 

The intergation over q in ( |S3|) can be done analytically, the one over g only 
numerically. However, under the simplifying assumption of a non-relativistic muon, 
i.e., M>p^>t;Cl^>g Cg, we can use the quasistatic approximation for the 
spectral functions 

cut, \ 3m 2 g 2 3m 2 g 1 

P? (Qo,q) = -^^\Dl(q ,q)\ 2 - 



2g 1 2g (g 2 + 3m 2 ) 2 ' 

cut, \ 3m 2 o;(g 2 - g 2 ) 3m 2 g g 

Pt (9b,?) = — — l^r(g ,g)| ~ 



4g 3 1 '^ u '*" 4 g 6 + (37rm2g /4) 2 ' 

Here we can easily identify the Debye screening, i.e. the screening of static electric 
fields associated with the Debye mass m 2 D = 3m 2 . Also we see that there is no static 
(?o = 0) magnetic screening. Using this approximation all integrations can be done 
exactly yielding 

e 2 T 



T L (v < 1) 



4tt ' 

e 2 T r dq 



Mv «d = e — v J 



? 
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As expected the longitudinal part, corresponding to the exchange of an electric pho- 
ton, is finite due to Debye screening, whereas the transverse part (exchange of mag- 
netic photons) is still IR divergent. Since q is restricted to soft momenta we adopt an 
UV-cutoff in the above expression for Ft of the order m 7 ~ eT. (The UV-behaviour 
1/q instead of 1/q 3 in Ft is a consequence of the quasistatic approximation.) The 
IR-cutoff cannot be calculated within the HTL method. However, one can show that 
the damping itself provides an IR-cutoff of order e 2 T. Then we obtain within the 
leading log approximation the final result 

e 2 T 1 

r T (v<l) = v\n-. (64) 

2n e 

Again let us make a few remarks here: 

1. We have seen that in the case of the damping rate the HTL method reduces the 
quadratic IR divergence, found in naive perurtbation theory, to a logarithmic one. 

2. Surprisingly T ~ e 2 and not of order e 4 as expected from naive perturbation 
theory. In other words, the HTL resummation leads to a lower order result (anomalous 
large damping) due to the strong IR sensitivty of the damping rate. 

3. In the case of a heavy quark in the QGP we simply have to replace e 2 in (|64]) 
by V/3. 

4. In QCD there is an alternative IR-cutoff, namely a magnetic screening mass of 
the order m magn = g 2 T, as suggested by non-perturbative arguments such as lattice 
QCD. 

4.2 Other quantities 

There is a number of other quantities of the QGP which have been computed using 
the HTL method. Here we want to give just a qualitative overview. For details we 
refer the reader to Refs.[2,3]. 

4.2.1 Energy loss, thermalization times, viscosity 

The energy loss of energetic quarks and gluons (E ^> T) in a QGP is related to jet 
quenching, which might serve as a signature for the QGP formation. Energetic quarks 
from initial hard collisions with a large transverse momentum have to propagate 
through the fireball. Depending on the phase the energy loss in the fireball might be 
different leading to a different energy distribution of jets. The energy loss per unit 
length is defined as the average of the energy transfer per collision AE divided by 




Fig. 30 
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the mean free path A = v/T. Hence it can be written as 

dE 



= V' 

v J 



. dTAE 

dx v J 

Due to the factor AE this expression is only logarithmically IR divergent in naive 
perturbation theory, i.e. 

dE r dq 
dx J q 

The energy loss can be calculated to leading order by introducing a separation scale 
gT <C q* <C T (g <C 1!). For q > q* it is sufficient to use a bare gluon propagator and 
we find from Fig. 29 



For q < q* we need a HTL gluon propagator as in Fig. 30 leading to 



' dE\ ri* a 

In 



rq 

/ dq... 
Jo 



*2 



\ dx I ,, jo ' 

\ / soft g 

Summing up the soft and the hard contributions we find 

dE _ fdE\ fdE\ ^ET 

dx \ dx )soit V^/hard 11 m 9 

where the arbitrary separation scale q* has dropped out. 

The complete calculation yields the generalization of the famous Bethe-Bloch 
formula to the case of energetic quarks in a QGP: 

dE 16vr „ - , 9E 
at T In ■ 



dx 9 s 16na s T 

shown in Fig.31 for an energetic charm and bottom quark. 

Note that dE/dx ~ a 2 s as in naive perturbation theory, because there is only a 
logarithmic IR singularity in naive perturbation theory in contrast to Y. Besides this 
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collisional energy loss due to elastic scattering, there is also a radiative energy loss 
due to gluon bremsstrahlung which turns out to be dominating in realistic situations. 

The thermalization time and the viscosity of the QGP are also closely related to 
the damping rate and can be calculated in a similar way as the energy loss, yielding 
for example for the thermalization or momentum relaxation time 



Since the energy in the logarithm - compare with the Bethe-Bloch formula - above is 
now replaced by a thermal energy the extrapolation from a s C 1 to a s > 0.2 breaks 
down. 

4.2.2 Photon and dilepton production 

The thermal emission of real and virtual photons from the QGP has been proposed as 
another promising signature for the QGP formation in relativistic heavy ion collisions. 
Let us first consider hard real photons with an energy E ^> T. To lowest order naive 
perturbation theory these photons are produced by the diagrams of Fig. 32. The 
photon production rate follows from the imaginary part of the two-loop photon self 
energy, which is related to the diagrams above again by cutting rules. In the case of 
a bare quark (massless) propagator we encounter a logarithmic IR singularity. Using 
the HTL method, we have to consider to lowest order the diagrams of Fig. 33, where 
one HTL quark propagator appears. Due to E 3> T it is sufficient to consider only 
one effective quark propagator. 

After some more or less tedious calculations one finds for the photon production 
rate 



which is shown in Fig. 34. 

Soft dileptons, i.e., lepton pairs from the decay of soft virtual photons (E, p ~ gT), 
can also be treated within the HTL resummation technique. Since the external photon 
momentum of the photon self energy is soft now, two effective quark propagators 
and also effective quark-photon vertices are necessary as in Fig. 35. According to 



r 



-i 





(65) 




Fig. 33 
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cutting rules this diagrams contains as physical process the annihilation of collective 
quarks and antiquarks. The non-trivial dispersion relation of the in-medium quarks 
in the HTL approximation leads to so-called Van Hove singularities in the dilepton 
production rate as in Fig.36, where a quantity proportional to Imll^^ is depicted. 

Recently it has been shown that there are additional contributions to the same 
order from higher order diagrams (e.g. bremsstrahlung) if the invariant photon mass 
of is of order g 2 T. Consequently the photon and ultrasoft dilepton production rates 
are non-perturbative and cannot be calculated to leading order even using the HTL 
improved perturbation theory. 

Let us shortly summarize the achievements of the HTL resummation technique: 

1. Using HTL resummed propagators and vertices for soft momenta, contributions 
of same (or even lower) order at T > as using bare propagators and vertices can be 
computed. 

2. At the same time medium effects (Landau damping, Debye screening) are 
included leading, e.g. to an improved IR-behaviour. 

3. Most important it provides gauge independent results for physical quantities 
in contrast to naive perturbation theory. 

However, there are still unsolved problems: 

1. The assumption g <C 1 is in contrast to realistic situations. In a QGP, expected 
to be produced in relativistic heavy ion collisions, we expect a s = 47i~/g 2 ~ 0.2 — 0.5, 
i.e. g = 1.5 — 2.5. The extrapolation of HTL results to such large couplings is 




Fig. 35 
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questionable and in the case of thermal energies E ~ T impossible (see e.g. the 
thermalization time). 

2. There is no static magnetic screening in the HTL approximation. Hence IR 
singularities from the magnetic sector are not removed completely (see e.g. damping 
rates). 

3. Some quantities cannot be calculated to leading order within HTL improved 
perturbation theory due to a strong IR sensitivity. For these quantities a non- 
perturbative resummation beyond HTL seems to be required (see e.g. the photon 
production rate). 

4. In relativistic heavy ion collisions the QGP fireball will be out of equilibrium at 
least in the early stage. Therefore an extension of the HTL method to non-equilibrium 
situations is desirable. 
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5 Non-Equilibrium 



If a partonic fireball is created in relativistic heavy ion collisions, it will be out of 
equilibrium at the beginning. Only by secondary collisions a thermal and chemical 
equilibrated QGP might be reached. Therefore we have to study non-equilibrium 
situations to describe the early stages of the fireball. 

Transport models, based on the Boltzmann equation for quarks and gluons, predict 
a fast thermalization (Thermal < 1 fm/c), but a slow (or no) chemical equilibration. 
Hence we have to consider non-equilibrium distributions for quarks and gluons. A 
simple ansatz for theses distributions in thermal equilibrium, but out of chemical 
equilibrium is the following: fp t B = A(i) riF,B(ko)- Here we simply multiply the equi- 
librium distributions by a time-dependent factor ( "fugacity" ) , A(t), which describes 
the deviation from chemical equilibrium. At the beginning there will be less quarks 
and gluons than in equilibrium, i.e. < A(t) < 1. 

Here we will discuss the generalization of the HTL method for quasistatic non- 
equilibrium situations, where the equilibration is slow compared to the processes 
under consideration^. We are not aiming at a description of the equilibration process, 
which is intrinsically non-perturbative. The starting point for this program is the 
RTF, since the ITF applies only to equilibrium. The basic idea is simply to replace 
the equilibrium distribution functions in (|25|) and fl32|) by non-equilibrium ones. Then 
it is easy to see that ([39]) to (|41|), i.e. the HTL results for the retarded and advanced 



polarisation tensor, still hold with 

m 7 -> m 7 = — / dkk f F (k) F -^X F (66) 

07T J o y 

where we assumed a locally isotropic momentum distribution fp = f F (k, x) depend- 
ing on k and the space-time coordinate x. Note that in deriving fl39|) to (41) no 



assumption about the existence of a temperature has been made. The HTL method 
relies only on the separation of scales. Instead of the temperature T one might use 
the average momentum (k). 

Instead of (^) (symmetric HTL polarization tensor) we get 



4?7> 2 r°° 

u L s (p) = ——e(p 2 -pl)J o dkk 2 f F (k)[i-f F (k)] 



with 



, Imn^(P) ,09 

= 2iA if pi < p 2 67 

Po 

^dkk 2 f F {k) [l-f F (k)] 

f °°dkkf F (k) • 1 J 

In equilibrium A reduces to 2T. Eq. ( p3|) (retarded and advanced HTL photon 

propagator) holds also using the replacement m 7 — > m 7 . 

4 The ideas presented in this section are based on Ref . [4] . 
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The most interesting quantity is the symmetric HTL propagator. We cannot use 
(|35| ) here anymore as we did for deriving (|54|). Rather we have to go back to the 
Dyson-Schwinger equation for D*g (see problem #9). It reads 

D* L S = D* L n + D% = + £ D^D% + D L 22 + £ D%$D% 

i,j=l i,j=l 



Using ( |29| ) and the inverse relations of (^) we find after some manipulations 

D* L S = D L S + D L R Tl L R D* L F + D L s Tl L A D* L A + D L R U L S D* L A . (69) 
It can be shown that this equation is solved by 

D* L S (P) = [1 + 2/b(po)] sgn(po) [I>*S(P) - D* L A (P)} 

+{n|(P) - [i + 2/ B (pb)] sgn(po) [n£(p) - ni(P)]} p>*£(p) d* l a (p). (70) 

Here we encounter a term ~ D* r{P)D* a{P)- In the case of bare propagators such 
a term produces a pinch singularity, i.e. terms containing S 2 (P). In equilibrium, where 
( Plf ) holds, the term in the curly brackets vanishes, {...} = 0. This is a consequence 
of detailed balance. Hence there is no pinch singularity in equilibrium. 

Now we study the dangerous pinch term, i.e. the second term of ([TfJ). The 
following equations hold in general: 

D*R, A = (p 2 -ReTl L R TiImTl L R )- 1 (71) 

and 

n£-ni = 2zimni (72) 

From (|7T|) we obtain 

* A (p 2 -Ren£) 2 + (Imn£) 2 " 2ilmn^ ' 1 ' 

if huH^ 7^ 0, i.e. for Pq < p 2 within the HTL approximation. Inserting (^) and ([73]) 
in (|70D yields 

s 2ilmU ] R K R AJ K J 

Hence there is no pinch singularity if Imn|j ^ 0. 

The physical interpretation is the following: in the case Imn|j ^ the particles 
have a finite width (damping). Consequently the 5 function is replaced by a func- 
tion p(P) of the Breit-Wigner form, i.e. a smeared 5-function. Using (|67|) and the 
definition of the spectral function (see section 3.2.2) one finds 

D* L S (P) = -2m - p L {P) (75) 
Po 

with Pl{P) follows from Pl{P) using the replacement m 7 — > m 7 . In equilibrium, 
where we have A — > 2T, (^) reduces to (|54|) in the HTL approximation p < T, for 
which [1 + 2n B (p )} sgn(p ) -> 2T/p . 
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Fig. 37 



As an application we discuss again the damping rate of a muon in a QED plasma. 
The leading order contribution in naive perturbation theory suffers from a pinch 
singularity in non-equilibrium since there are two photon propagators of the same 
momentum (see Fig. 37). In the HTL approximation, however, the 2-loop diagram of 
Fig. 38 is replaced by a 1-loop diagram containing an effective gluon propagator, as 
in section 4.1. The photon momentum is restricted to the Q 2 < (see (|63|)), where 
Imll^ ^ 0. Hence the HTL method also removes the pinch singularity found in naive 
perturbation theory. Eq. ( |B"3"D reads now 

pl(qo, q) + [v 2 - ^ priio, q) 

(76) 

This result does not depend on m 7 since it drops out after performing the integrations 
in©- 

Summarizing, we have shown, how the HTL technique can be generalized to qua- 
sistatic non-equilibrium situations. The symmetric HTL propagator contains no pinch 
singularity below the light cone due to Landau damping. In general, we can say that 
pinch singularities come from using bare propagators in perturbation theory. Re- 
summed propagators for quasiparticles with a finite width (Breit-Wigner) do not 
lead to pinch singularities. In equilibrium pinch singularities are absent also for non- 
interacting particles due to detailed balance. 



r neq (E) = - — / dqq / 

4ttv Jo J-vq q 

= r eq (£). 




Fig. 38 
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Problems 



1. Calculate the particle and energy density of an ideal quark-gluon plasma con- 
taining up-, down-quarks, and gluons. Assume that the quarks are massless 
and that there are as many antiquarks as quarks. Use the result to estimate the 
critical particle and energy density for the deconfmement transition at T c = 170 
MeV. 



Hint: Use 



(n-1)! Jo e x -l (1 - 2 1 - n )(n - 1)! Jo e x + 1 
with C(3) = 1.202 and C(4) = tt 4 /90. 



2. Show that 



1 f n 

I = T £ f(k = 2mnT) = — / 



1 f r° JL f(ko)+f(-k ) 



dko 



•JOG 

ioo+e 



+ 



/ioo+e ~i 
dfco[/(fco) + /(-fco)]n B (A;o)[, 
-ioo+e J 



if /(/co) has no poles on the imaginary axis. 
Hint: Show first 

I = ^-<f c dk f(k )coth(f3k /2), 
where C are the contours around the poles of coth(/?fc /2). 

3. Show that the propagator in the Keldysh representation 

( Aa \ 



A K -- 

follows from the RTF propagator 



A 5 



/A n A 12 
\ A 21 A 22 



by applying an orthogonal transformation 
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4. Show that the symmetric propagator can be written as 

A S (K) = [1 + 2n B (k )} sgn(A;o) [A R (K) - A A (K)}. 



5. Derive the symmetric longitudinal photon self energy in the HTL approximation 
and show that it fulfils the relation 

U L S (K) = [1 + 2n B (k )] sgn(k ) [U L R (K) - U L A (K)]. 

Hint: Use 

n F (k)[l - n F (k)\ = -T—n F (k). 



6. Show that the Dyson-Schwinger equation for the full photon propagator in 
Coulomb gauge leads to 



7. Calculate the photon dispersion relation in the HTL approximation for small 
k <C uj up to order k 2 . 

Hint: In i±£ = 2(x + *f + ^ + ...)• 

8. Show that the tadpole diagram in massless 4 -theory containing the effective 
propagator with mass M = gT (sum of daisy diagrams) is given for small g< 1 
by 

U* = g 2 T 2 \l-- g + 0(g 2 ) 

Hint: Show that the UV divergent part is of 0(g 4 ) after regularizing it by 
subtracting the 1-loop tadpole at T = and introduce a separation scale gT 
k* < T for evaluating the momentum integral in the small g limit. 

9. Show that the Dyson-Schwinger equation for the full symmetric propagator 
reads 

A* S = A S + A R U R A* S + A S U A A* A + A R U S A* A 
and that it is solved by the ansatz 
A S (K) = [l + 2n B (k )]sgn(k )[A R (K)-A A (K)] 

+ {U S (K) - [1 + 2n B (k )]sgn(k )[Il R (K) - U A (K)]}A R (K) A* A (K). 

Hint: Use A^n^A^ = A* R A - A R:A . 
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